Ultracold Bose Gases in ID Disorder: From Lifshits Glass to Bose-Einstein Condensate 
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We study an ultracold Bose gas in the presence of ID disorder for repulsive inter-atomic interactions varying 
from zero to the Thomas-Fermi regime. We show that for weak interactions the Bose gas populates a finite 
number of localized single-particle Lifshits states, while for strong interactions a delocalized disordered Bose- 
Einstein condensate is formed. We discuss the schematic quantum-state diagram and derive the equations of 
state for various regimes. 
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Disorder is present in nearly all condensed-matter systems 
due to unavoidable defects of the sustaining media. It is 
known not only to impair quantum flows but also to lead to 
spectacular effects such as Anderson localization ([H 0. Ell- 
in contrast to condensed-matter systems, ultracold atomic 
gases can be realized in the presence of controlled disor- 
der or quasi-disorder |4|], opening possibilities for investiga- 
tions of localization effects (for review see 
Ref. Gob. The first experimental studies of localization 
in disordered interac ting Bose gases have been reported in 
Refs. ifTTl [l2l [T3l fl4l |T5r[T6tl ■ 

One of the most fundamental issues in this respect concerns 
the interplay between localization and interactions in many- 
body quantum systems at zero temperature. Without interac- 
tions, a quantum gas in a random potential populates local- 
ized states [1], either a single state (in the case of bosons), or 
many (fermions). Weak repulsive interactions lead to der- 
ealization but strong interactions in confined geometries lead 
to Wigner-Mott-like localization 111 711 . Surprisingly, even for 
weakly interacting Bose gases, where the mean-field Hartree- 
Fock-Gross-Pitaevskii-Bogolyubov-de Gennes (HFGPBdG) 
description is expected to be valid, there exists no clear picture 
of the localization-delocalization scenario. Numerical calcu- 
lations using the Gross-Pitaevskii equation (GPE) suggest that 
the Bose gas wave-function at low densities is a superposi- 
tion of localized states lfl5ll . It is thus natural to seek the 
true ground state in the form of generalized HFGPBdG states, 
for which the Bose gas populates various low-energy single- 
atom states. In the presence of disorder, they correspond to 
so-called Lifshits states (LS) [18]. 

In this Letter we consider a d-dimensional (dD) Bose gas at 
zero temperature with repulsive interactions, and placed in a 
ID random potential with arbitrary amplitude and correlation 
length. We show that generalized HFGPBdG states indeed 
provide a very good description of the many-body ground 
state for interactions varying from zero to the Thomas-Fermi 
(TF) regime fl9ll . We stress that the solution we find is differ- 
ent from that of non-interacting fermions which at zero tem- 
perature form a Fermi glass and occupy a large number of 
localized single particle levels lioll . In contrast, many bosons 
may occupy the same level and thus populate only a finite 



number of LSs forming what we call a Lifshits glass. In 
the following, we discuss the quantum states of the system 
as a function of the strength of interactions and the ampli- 
tude and correlation length of the random potential, and we 
draw the schematic quantum-state diagram (see Fig. [TJ. In 
the limit of weak interactions, the Bose gas is in the Lifshits 
glass state, whereas for stronger interactions the gas forms 
a (possibly smoothed) delocalized disordered Bose-Einstein 
condensate (BEC) 12011 . Our theoretical treatment provides us 
with a novel, physically clear, picture of disordered, weakly- 
interacting, ultracold Bose gases. This is the main result of 
this work. In addition, we derive analytical formulae for the 
boundaries (corresponding to crossovers) in the quantum-state 
diagram and for the equations of state in the various regimes. 
We illustrate our results using a speckle random potential 12 ill . 

Consider a dD ultracold Bose gas with weak repulsive inter- 
actions, i.e. such as n 1 ~ 2 / d -C h 2 /mg, where m is the atomic 
mass, n the density and g the dD coupling constant. The gas 
is assumed to be axially confined to a box of length 2L in the 
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Figure 1: (color online) Schematic quantum-state diagram of an in- 
teracting ultracold Bose gas in ID disorder. The dashed lines rep- 
resent the boundaries (corresponding to crossovers) which are con- 
trolled by the parameter q r = H 2 /2mo^ Vr (fixed in the figure, see 
text), where Vr and cr R are the amplitude and correlation length of 
the random potential. The hatched part corresponds to a forbidden 

ZOne (jU < Vniin). 
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presence of a ID random potential V(z), and trapped radially 
in a 2D harmonic trap with frequency ui±. We assume that 
the random potential is bounded below [V m i„ = min(V)] and 
we use the scaling form V(z) = Vrv(z/<t r ), where v(u) is a 
random function with both typical amplitude and correlation 
length equal to unity I22I1 . 

For illustration, we will consider a ID speckle potential 
lElll similar to that used in Refs. (H El Oj, Q . In brief, 
v(u) is random with the probability distribution V{v) = Q(v+ 
1) exp[— (v + 1)], where is the Heaviside function. Thus 
v is bounded below by v m \ n = —1 and we have (v) = and 
(y 2 ) = 1. In addition, for a square aperture, the correlation 

]2(u-u') 



function reads (v(u)v(u')) =sinc 

Below, we discuss the quantum states of the Bose gas, 
which are determined by the interplay of interactions and dis- 
order. 

BEC regime - For strong repulsive interactions the Bose 
gas is delocalized and forms a BEC lfl2l I20I1 (possibly quasi- 
BEC in ID or elongated geometries 02311 '). The density profile 
is then governed by the GPE, 

(1 = -n 2 V 2 (V?i)/2mVn + mcj 2 _p 2 /2 + V(z) + gn(r), (1) 

where p is the radial coordinate and p the chemical poten- 
tial. This regime has been studied in the purely ID case in 
Ref. I20I1 . Here, we focus on the case of a shallow radial trap 
(hui±-^ip) such that the radial profile is a TF inverted parabola. 
Proceeding as in Ref. 12011. we find that the BEC density has a 
generalized TF profile 12411 : 



y/n(p, z) ~ y/fi(p)/g [1 - V(p, z)/2p(p)} (2) 

where p{p)=p\l— (p/R±) 2 ] is the local chemical potential, 
R± = i/2/i/mw 2 _ is the radial TF half-size and V_[ 
Jdz'G(p, z')V(z — z') is a smoothed potential 12C 



TtiP) ° XP IT ' and ^ = ^ V 2m M(^) be- 
ing the local healing length. For £(/9)<Ccr R , i.e. for 



G( P ,Z): 



fi(p) > h 2 /2ma R , 



(3) 



we have V(p, z) ~ V(z), and the BEC density follows the 
modulations of the random potential in the TF regime. For 
£(/9)><t r , the kinetic energy cannot be neglected and competes 
with the disorder and the interactions. The random potential 
is therefore smoothed lHofc: AV(p)<AV where AV (AV(p)) 
is the standard deviation of the (smoothed) random potential. 
The solution (ffji corresponds to a delocalized disordered BEC. 

The perturbative approach is valid when p(p) 3> 
AV{p). From the expression for V, we write AV(p) = 
Vr-\/S (ct r /^(p)). For the speckle potential, we can approx- 
imate the correlation function to exp (— z 2 /2cr 2 ) and we 
find 01 
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Figure 2: (color online) a) Cumulative density of states of single par- 
ticles in a speckle potential with a R = 2 x 10 _3 L and Vr = lO 4 ^), 
where Eq = h 2 /2rnL 2 (T4ii n = — Vr). Inset: Participation length 
I25I1 . b) Low-energy Lif shits eigenstates. For the considered realiza- 
tion of disorder, eo — —5 x W 3 Eg. 



with cr R =cr R /£(p). In the center, /.e. p = 0orin ID, the validity 
condition of the BEC regime thus reduces to 



H » V R VS°(a R /C) with C = C(0). 



(5) 



S (5 R ) =a 2 + (l-25 2 )a R e 3 



d6»e" 



(4) 



If condition (0 is not fulfilled, the Bose gas will form a 
fragmented BEC. The latter is a compressible insulator and 
thus can be identified with a Bose glass IU7I1 . 

Non- interacting regime - In the opposite situation, for van- 
ishing interactions, the problem is separable and the radial 
wave-function is the ground state of the radial harmonic os- 
cillator. We are thus left with the eigenproblem of the single- 
particle ID Hamiltonian h=— h 2 d 2 /2m + V(z). In the pres- 
ence of disorder, the eigenstates x„ are all localized J2] and 
are characterized by II 1 811 (i) a finite localization length, (ii) 
a dense pure point density of states T>2L, and (iii) a small 
participation length P v = 1/ Jdz \xv(z)f iSl- If V(z) is 
bounded below, so is the spectrum and the low energy states 
belong to the so-called Lifshits tail, which is characterized by 
a stretched exponential cumulative density of states (cDOS), 

AA 2L (e)=/ e de'P 2L (e')-cxp(- Cv ^^), in ID & 

Numerical results for the speckle potential are shown in 
Fig. |2 As expected the cDOS shows a stretched exponen- 
tial form, the lowest LSs are spatially localized, and P(e) in- 
creases with energy indicating a weaker localization. How- 
ever, P(e) is almost constant at low energy. Note also that the 
lowest LSs hardly overlap if their extension is much smaller 
than the system size. 



3 



Lifshits regime - We turn now to the regime of finite but 
weak interactions, where the chemical potential p lies in the 
Lifshits tail of the spectrum. Owing to the fact that the lowest 
single-particle LSs hardly overlap, it is convenient to work in 
the basis of the LSs, {xv, v 6 N}. These can be regarded as 
trapping micro-sites populated with N v bosons in the quantum 
state 4>v{p)Xv( z ) where the longitudinal motion is frozen to 
Xu an d 4>y accounts for the radial extension in the micro-site 
v. Therefore, the many-body wave-function is the Fock state 



l*> = n^O-^^lvac) 



(6) 



where b\ is the creation operator in the state <fiu(p)xv(z) 12711 . 
Each <j) v can be a transverse 2D BEC for N V ^>1. However, the 
quantum state (O does not correspond generally to a single 3D 
BEC since it does not reduce to (AH) -1 / 2 (bl) N |vac) . Rather, 
the Bose gas splits into several fragments whose longitudinal 
shapes are those of the LSs, %„, and are hardly affected by the 
interactions. 

The mean-field energy associated to the state © reads 

n 2 r 



(7) 



where U v =g jAz\xv{z)\ i =gP~ 1 is the local interaction en- 
ergy in the LS Xv Minimizing E[^>] for a fixed number of 
atoms (E[if?] — fi^2 N v — > min), we find the equation 

(p-e v )(f> v = [-h 2 ^ 2 j2m + mu\p 2 12 + N V \J V \$ V \ 2 \ <?V 

(8) 

Solving the 2D GPE ((HJ for each micro-site v, one finds the 
atom numbers N v and the wave-functions <f> v . As p increases, 
1 4> v \ 2 will turn continuously from a Gaussian (for huj±^>p) 
into an inverted parabola (for huj±<^p). 

To discuss the validity condition of the Lifshits regime, let 
us call v mdK the index of the highest LS such that all lower LSs 
hardly overlap. The Lifshits description requires the chemical 
potential p to be small enough so that the number of populated 
LSs is smaller than z/ max , i.e. if 



M 2 l{p) < v n 



(9) 



If condition (O is not fulfilled, several populated LSs will 
overlap and the Bose gas will start to form a fragmented BEC. 
Each fragment will be a superposition of LSs, and its shape 
will be modified by the interactions. 

Although both Mil an d ^ max ma y have complex dependen- 
cies versus Vr, <j r and the model of disorder, general proper- 
ties can be obtained using scaling arguments. We rewrite the 
single-particle problem as 



(£„/Vr)vj„(u) = ~a R d 2 (p v (u) + v(u)tp v (u), 



(10) 



where u~z/a R , tp v (u) = ^/a^x^( z ) an d a R — fi 2 /2ma 2 VR. 
Thus, all characteristics of the spectrum depend only on the 



parameter a R after renormalization of energies and lengths. 
Scaling arguments show that in the Lifshits tail 

Af 2L (e) = (L/a s )((a R ,e/V R ) and v'™ = (L/a R )r,(a R ), 

(ID 

where £ and rj are w-dependent functions. Finally, inserting 
these expressions into Eq. (O and solving formally, we obtain 
the validity condition of the Lifshits regime: 



P < V R F{a R ), 



(12) 



where F is the solution of ((a R , i^(a R ))=7y(a R ), which can be 
computed numerically, for example. 

We are now able to draw the schematic quantum-state dia- 
gram of the zero-temperature Bose gas as a function of p and 
Vr (see Fig.HJ. From the discussion above, it is clearly fruit- 
ful to fix the parameter a R while varying Vr. The boundaries 
between the various regimes (Lifshits, fragmented BEC, BEC 
and smoothed BEC) result from the competition between the 
interactions and the disorder and are given by Eqs. (1315112b . 
We stress that they are crossovers rather than phase transi- 
tions. Interestingly, all these boundaries are straight lines with 
slopes depending on the parameter a E . This is clear from 
Eq. ( fT2b for the boundary between the Lifshits and the frag- 
mented regimes. In addition, since VR~(h 2 /2ma 2 )/a R , the 
non-smoothing condition (O reduces to p ^> C( r Vr. Finally, 
since p = h 2 /2to£ 2 and thus ct r /£ = — h^\J p/Vr, the non- 
fragmented BEC condition <(5j also corresponds to a straight 
line with a slope depending on a R in Fig.Q] 

To finish with, we derive the equations of state of the Bose 
gas in the identified quantum states. It is important to relate 
the chemical potential p which governs the crossovers be- 
tween the various regimes to the mean atomic density n = 
N/2L and the coupling constant g. Both can be controlled in 
experiments with ultracold atoms. 

Tight radial confinement - For p'—e u <tifruJ±_ where p'=p— 
hoj±, the radial wave-functions are frozen to zero-point oscil- 
lations, 4>v{p) = cxp(—p 2 /2l 2 J/y/Trl± with l± = yjh/muj± 
the width of the radial oscillator. 



In the BEC regime, p! 3> AV^, we find from Eq. (0, 



p! = ng. 



(13) 



In the Lifshits regime, we find 



N v =[p' — e„]/E/„ for p! > e v and N U = Q otherwise, (14) 

by inserting the above expression for 0„ (p) into Eq. ([8]). Turn- 
ing to a continuous formulation and using the normalization 
condition, N = JdeT>2L{^)N(e), we deduce the equation of 
state of the Bose gas in the Lifshits regime: 



Ng 



de V 2L (e) (y-e)P(e), 



(15) 



which relates the chemical potential p! to the coupling con- 
stant g. The relation is in general non-universal (i.e. it depends 
on the model of disorder, v). In the case of a speckle potential, 



4 



25000 
20000 
15000 
10000 
5000 


-5000 
-10000 



— Lifshits regime 

■■■ B EC regime 

• numerical meanfield 



* BEC 

regime 




10000 15000 
ng/E 



20000 25000 



formulae for the boundaries (crossovers) in the quantum-state 
diagram and shown that they are determined by the coupling 
constant. Since this coupling constant can be controlled in 
cold gases, future experiments should be able to explore the 
whole diagram. 
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Spanish MEC Programs, DUQUAG, QOIT and the European 
Union QUDEDIS and IP SCALA. 



Figure 3: (color online) Chemical potential of a Bose gas in a speckle 
potential with the same parameters as in Fig. [2] in the case of a tight 
radial confinement (// — e v <C fujf). The points are given by numer- 
ical calculations; the solid and dotted lines represent the analytical 
formulae Eqs. dl6t and J 1 3t > derived for the Lifshits and BEC regime, 
respectively. 



■A/2i,(e)=A(Q! R )(i/(j R ) exp (-c(a R ) / ^ e/V R +l) , and assum- 
ing that the participation length P v =cr R p u (a R ) is independent 
of the energy in the Lifshits tail, we find: 

rig ~ A(a R )c 2 (a R )p (a R )V R T (-2, c(a R )/yJ p' /V R + ij , 

(16) 

where T is the incomplete gamma function and A(a R ), c(a R ) 
and po(a R ) can determined numerically. 

Using a numerical minimization of the energy func- 
tional (|7|i in the Gross-Pitaevskii formulation, we compute the 
chemical potential of the Bose gas in a wide range of interac- 
tions. The result shown in Fig. [3] indicates a clear crossover 
from the Lifshits regime to the BEC regime as the interaction 
strength increases. The numerically obtained chemical poten- 
tial [i agrees with our analytical formulae in both Lifshits and 
BEC regimes. 

Shallow radial confinement - The equations of state can 
also be obtained in the case of shallow radial confinement 
(fi — e„ S> huj±). In the BEC regime, for £ <C <J R , we 
find fi ~ yjngmuj 2 ^ / tt— . In the Lifshits regime the 2D 
wave-functions 4> v {p) are in the TF regime, \4>»(p)\ 2 = 
§^ (l-p 2 /Rl), where R u = y/2([i.-eJ)/mu{ is the 2D- 
TF radius and N v = ■n{p — e u ) 2 /mw 2 _U u for fj, > e u (0 other- 
wise). Proceeding as in the ID case, we find: 

Ng = ^ I deV 2L (e)( f i-efP(e). (17) 

muj i J-oo 

Applying this formula to the relevant model of disorder allows 
us to compute the populations N v of the various LSs \v an d 
the corresponding radial extensions <j> v . 

In summary, we have presented a complete picture of the 
quantum states of an interacting Bose gas in the presence of 
ID disorder, including the novel description of the weakly 
interacting Lifshits glass state. We have provided analytical 
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